10 /tW 
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/define MIN_HI 1 
#define MAXJil 2001 
#def ine HI_INDEX 1 

#define MAX_IDX 4000 

Tissue-equivalent epoxy plastic 
C4 plastic explosive */ 
RDX sheet explosive */ 
Water Gel explosive */ 
40% dynamite stick */ 



/* #define 


TISSUE */ 


/* 


/define C4 




/* 


/* #define 


RDX */ 


/* 


/* /define 


WG */ 


/* 


/* #define 


DYN */ 


/* 



/* new way of determining low */ 

#define zl .0247 

#define z2 .01492 

#define z3 .265 

/define z4 112.6 

#define z5 25.198 

/define z6 .6218 

#define z7 .265 

/* define substance parameters */ 



#ifdef WG 
#define cl 9.732 
#define c2 6.108 
/define c5 1.218 
/define K0 .547 
/define KL .961 
/endif 



/ifdef RDX 
/define cl 9.732 
/define c2 6.108 
/define c5 1.218 
/define K0 .65 
/define KL .86 
/endif 



/ifdef C4 
/define cl 9.732 
/define c2 6.108 
/define c5 1.218 
/define K0 .6522 
/define KL .87 
/endif 

/ifdef DYN 
/define cl 570.46 
/define c2 4.352 
/define c5 .304 
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#define K0 .522 
#define KL .765 
#endif 



#ifdef TISSUE 
#define cl 3798 
/define c2 3.8837 
#define c5 0.993 
/define KO .655 
/define KL .825 
#endif 

double bh (double km) ; 
- double bh (double km) 



double Kref (double Hi, double Km, double kO) ; 
double Kref (double Hi, double Km, double kO) 

{ 

return ( ( (Hi+bh (Km) ) *kO*KL) / ( (bh (Km) *KL) + (Hi*ko) ) ) ; 

} 



double alpha (double km); 
double alpha (double km) 

{ 

retum((zl+(z2*km)-(z2*z3) )/(km*km) ) ; 

> 

double beta (double km) ; 
double beta (double km) 
{ 

return ( (z4+( (z6-km) * (z5/ (z6-z7) ) ) ) /km) ; 

} 

double newlow (double h, double km); 
double newlow (double h, double km) 

{ 

return (h*(l/ (km+ (alpha (km) * (h/ (h+beta(km) ) ) ) ) ) ) ; 



double find_Km (double hi, double Kair,double kref); 
double find_Km( double hi, double Kair, double kref) 

/* find the Km that approximates the desired Kref given high val,k0 */ 
int x,bitval; 

double lsbval , approx_kref ; 



{ 



return (cl*pow( (km+c5) ,c2) ) ; 



} 



} 
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lsbval = 0.8; 
bitval =0; 

for (x=0;x<8;x++) 

{ bitval=(bitval«l) 1 1; 

lsbval = lsbval/ (double) 2.0; 

approx_kref = (Kref (hi, ( (double) .l+( (double) bitval*lsbval) ) ,Kair) ) ; 
if (approx_kref < kref) 



} 

return ( ( (double) bitval*lsbval) +. 1) ; 

} 

double findKm_Low( double hi, double low); 
double findKm_Low (double hi, double low) 
{ 

/* find the Km that approximates the desired Low given high val,k0 */ 

int x, bitval; 

double lsbval , approx_low ; 

lsbval = 0.8; 
bitval = 0; 

for (x=0 ; x<8 ; x++) 



bitval=(bitval«l) |l; 
lsbval = lsbval/ (double) 2.0; 

approx_low - (Low (hi, ( (double) . 1+ ( (double) bitval* lsbval) ) ) ) ; 

if (approx_low < low) 

bitval=bitval&(0xfe) ; 



return ( ( (double) bitval*lsbval)+.l) ; 



/* create the histogram */ 

for (hint = MIN_HI; hint < MAX_HI; hint += HI_INDEX) 
{ 

h = (double)hint; /* Get hi double valije . */ 



bitval=bitval& (Oxf e) 



{ 



} 



/* Set up the header values and the KIdx */ 
Hdr[HI_VALUE] = hint; 
KIdx = 0; 
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/* Get the hi and lo kref */ 
hikref = Kref(h, 0.29, kO) ; 
lo_kref = Kref(h, 0.8, kO) ; 
k=lo kref; 



lastl = -100.0; 
diffl = 1000.0; 
while (k<hi_kref) 
{ 

km=f ind_Km (h , kO , k) ; 
kr=Kref (h,km,kO) ; 
l=Low(h,km) ; 

if (((l-lastl)<diffl)&&(km>.29)) 
diffl = 1 - lastl; 

lastl = 1; 

if (h>800.0) 
{ 

k=k*1.04; 
} else 

k=1.01*k; /* 1% bins */ 



/* do it again, but use diffl to find values */ 

k=lo kref; 

km=fTnd_Km(h,k0,k) ; 

l=Low(h,km) ; 

findl=(int)l; 



/* adjust diffl to a power of 2 */ 
tdiffl=0; 

while ((1 « (tdiffl+1)) <= (int) diffl) 
tdiffl++; 

km=f indKm_Low(h, (double) findl) ; 
k=Kref (h,km,k0) ; 

/* Save the minimum low and the scale factor */ 
Hdr [MIN_L0] = findl; 
Hdr [ LO_SCALE ] = tdiffl; 

while (k < hi_kref ) 
{ 

km=f indKm_Low (h , ( double ) f indl ) ; 
k=Kref (h,km,k0) ; 

/* Save the necessary information into the values */ 
KrefTab[KIdx] = (float) k; 
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KIdx++; 

/* increment low */ 

f indl += (1 « tdif f 1) ; 

/* increment bin count */ 
bincnt+=l; 



/* Now we have the table, write out the header then the table */ 
Hdr[MAX_LO] = findl; 

bwritten = write (fhndl, (char *)Hdr, sizeof (int) *4) ; 
if (bwritten != (sizeof (int) * 4)) 

* pr int f ("Error writing f ile\n") ; 
return (1) ; 

} 

/* Now write out the kref vector */ 

bwritten = write (fhndl, (char *)KrefTab, sizeof (float) *KIdx) ; 

if (bwritten != (sizeof (float) *KIdx) ) 

{ 

printf ("Error writing f ile\n") ; 
return (1) ; 

} 



/* output bin count */ 

printf ("Total Kref bin count :%ld\n",bincnt) ; 



} 



} 



/* 



Detection algorithm for above histogram 



*/ 



/* 



* Function: 

* DoBox 



* Descrition: 



Process the box. 



* Usage: 



DoBox (x, y) 



* Inputs : 



x - int : the x coordinate of the candidate pixel 
y - int : the y coordinate of the candidate pixel 




t 



t 
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* Outputs: 

* None 



*/ 



static void DoBox (int x, int y) 
{ 

int tx, ty; 

double diffH, diffL, diffK; 
double kref lo / krefhi ,krefavg; 
/* int tmp; */ 
double mindiff; 

Pixel *midpxl = &ScanLine[y] [x] ; 
Pixel *pxl; 

/* Average the values for this pixel */ 
AveragePixel (x, y) ; 

/* See if we need to do this pixel */ 
if (midpxl->avghia > 2000.0) 



* Calculate the min difference value (this is calculated by using 

* twice the expected noise as the difference value) . 



mindiff = (10000. 0/ (100. 0+midpxl->avghia) ) ; 

/* Now loop through the pixels doing the box */ 
for (ty = y - BORDER; ty <= (y + BORDER) ; ty++) 
{ 

/* Get the pixel */ 

pxl = &ScanLine[ty] [x - BORDER]; 

/* Loop through the x */ 

for (tx = x - BORDER; tx <= (x + BORDER) ; tx++, pxl++) 
{ 



/* See if we need to look at this pixel (edges are no-nos) */ 
if (pxl->sobel) 
continue ; 

/* Average this sucker */ 
AveragePixel (tx # ty) ; 

/* Now difference the Hi AIRS */ 
diffH = midpxl->avghia - pxl->avghia; 

/* Now threshold it */ 
if (diffH < mindiff) 
continue; 

/* Now difference the Lo AIRS */ 
diffL = midpxl->avgloa - pxl->avgloa; 



return ; 



/* 



*/ 




* 



• 



ALGO.DOC 



Thursday, May 3, 1990 9:12 am 



Page 7 



} 

if 



/* Now threshold it */ 

if ((diffL < mindiff) || (diffL == 0.0)) 
continue; 

kref lo=LookupKref (pxl->avghia,pxl->avgloa) ; 

kref hi=LookupKref (midpxl->avghia,midpxl->avgloa) ; 

diffK = diffH/diffL; 

/* Key lookup algorithm 
* Histogram generation algorithm has been fit to this ratio 

*/ 

kref avg= (kref lo* . 8) + ( . 2*kref hi) ; 

/* See if we need to histogram this point */ 
if ((diffK < (krefavg+(MinThreshold) ) ) 

| | (diffK >(krefavg+(MaxThreshold) ) ) ) 

continue; 

midpxl->histval++ ; 

} 



(maxhit<midpxl->histval) 

maxhit==midpxl->histval ; 
midpxl->histval > fomThresh) 

f om += (midpxl->histval - fomThresh) ; 

( (midpxl->histval > 0) && (midpxl->histval <200) ) 
histpix [midpxl->histval ] ++ ; 



} 




